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Phase determination from mostly one-sided interferograms 
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Abstract 

We show how to detect and correct for non-linear phase shifts in a mainly one-sided 
interferogram of an emission-line source. We simultaneously detect and correct for 
an out-of-phase emission background from the spectrometer. The method requires 
two auxiliary spectra, one of a strong continuum source, and one of an emission-line 
source with little or no continuum. 
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1. Introduction 

Phase determination for an ideal interferometer amounts to finding the position of zero path 
difference (ZPD) for each spectral element in a recorded interferogram. In this case, the 
interferogram F(x') of a single unresolved spectral line centered at the wavelength A 0 = l/cr 0 
is given by F(x') = a cos 27roo(x' —x), where x is the position of ZPD and x' is the retardation 
of the retroreflector in the long arm of the interferometer. The instantaneous optical path 
difference (OPD) is given by x' — x. In the ideal case where an interferogram sample point 
falls precisely on the peak of the interferogram we have x = 0, but in general x ^ 0. The 
transform of a one-sided interferogram gives 

F(a) = \ [ L F(x') e iW dx' 
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where L is the length of the scan, a = 2n (a — ao)L, and we have dropped terms containing 
(<7 + cr 0 ). The transform of a mostly one-sided interferogram (where the integral limits are 
— L\ to +L 2 , Li <C L 2 , and a tapering function is applied between —Li and +L\ to avoid 
double-counting of the two-sided part of F) is similar, but differs in details which are not 
important here. In either case, once we have determined x the spectrum is recovered by 
taking the real part of F(o)e~ l2 ™ 0X . The transform of a two-sided interferogram (where the 
mirror scans from — L to L) for the same unresolved line is simply ae l2n,JoX sin a /2a. 

The phase of the transform F(a) in Eq. 1 is given by a slowly varying term (27ra 0 x), plus 
a rapidly varying term, tan -1 [(1 — cos a ) / sin a], which is only significant near a spectral line. 
We have found that, for our interferometer, the slowly varying term is not well modelled by 
the linear function 2i rax. In Sections 3-5 we present a method for detecting and correcting 
for a small non-linear phase shift. We simultaneously detect and correct for an out-of-phase 
emission background from the spectrometer. 

Eq. 1 illustrates a crucial element of our phase-recovery method, line symmetrizing, 
which is discussed in Section 4. The phase factor Ittgqx in Eq. 1 is not known a priori , but 
must be estimated. If the estimated phase equals 2-koqx, then the real part of Eq. 1 gives 
a symmetric line, from the sin a/ a term. If the estimated phase is 2n<7 0 x ± 7t/2, then the 
real part of Eq. 1 gives an antisymmetric line, from the (1 — cosa)/o; term. In practice, it 
is straightforward to select the optimum phase factor which symmetrizes the real part of a 
complex spectral line by multiplying the region around the line by a series of equispaced 
phase factors, and testing the real part of each result for symmetry. 

2. Instrument description 

The far-infrared spectrometer (FIRS-2) was designed and built at the Smithsonian As- 
trophysical Observatory for the purpose of measuring stratospheric emission spectra from 
balloon and aircraft platforms. The spectrometer produces mostly one-sided interferograms, 
which must be transformed and accurately phase-corrected before the spectra can be used 
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to retrieve constituent profiles. The FIRS-2 and data reduction procedure are described in 
detail elsewhere 1,2 ; our phase-recovery method is summarized in the latter 2 , but fully dis- 
cussed in the present paper. A brief description of the instrument as configured for balloon 
flights follows. 

The flight instrument consists of a limb-scanning telescope, the interferometer, detectors, 
and associated control electronics. The telescope is a small off-axis reflector which is pointed 
in elevation relative to a single-axis stabilized platform 3,4 . The interferometer uses hollow 
corner-cube retroreflectors and all pm uncoated Mylar beamsplitter, and scans over optical 
path differences from -L x to +L 2 , where L x = 1.2 cm and L 2 = 120 cm. We use two liquid 
helium cooled detectors, a gallium-doped germanium photoconductor for the 75-220 cm -1 
band and a copper-doped germanium photoconductor for the 330-700 cm -1 band. The field 
of view is defined by field stops near the detectors. 

During a balloon flight we record interferograms for a sequence of elevation angles, repeat- 
ing the sequence throughout the flight. We first observe an ambient-temperature blackbody 
source at an elevation angle of 90°, followed by a space view at 30°, and then we observe 
the limb at elevation angles of 0.0, -2.3, -3.0, -3.6, -4.2, and -4.6°. A complete sequence 
takes about 2700 s. We show small segments of sample flight spectra for 90, 30, and -4.17° 
in Fig. 1, recorded at a balloon altitude of 36 km. The spectra have been phase corrected 
using the method outlined below, but have not been calibrated in intensity or wavenumber. 
Note the presence of a substantial out-of-phase (negative) instrumental background in the 
30° spectrum. The phase recovery method discussed here is not affected by this background, 
which is later detected and corrected for. 

3. Model description 

We assume that the measured blackbody spectrum Mg 0 (a) can be modelled as 

M 90 (<t) = [B{a) + (2) 
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where B(a) is the product of the instrumental response function and a blackbody spectrum, 
E(a)e l<t> is the complex background spectrum 5 , £90 is the position of ZPD in the blackbody 
interferogram, and e(a) is the non-linear phase shift. 

Similarly, we model the measured space and limb spectra as 

M 30 (a) = [S 30 (a) + (3) 

M g {a) = [S*(<r) + (4) 

where S 30 and Sg are the phase-corrected space and limb spectra, and x 30 and xg are the 
positions of ZPD in the respective interferograms. 

In this notation, the A/go^o^c) functions are complex spectra obtained by applying 
a fast-Fourier transform to an observed and windowed interferogram, where the window 
function is typically a weight running linearly from 0 to 1 between —Li and +L\, and a 
constant weight of 1 between +Li and +L 2 . Also in this notation, the functions on the right 
sides of Eqns. 2-4 ( B(cr ), £ 30 , ^(cr), E(cr), <t>(cr), e(cr), and £ 90 , 30 , 0 ) are all to be extracted from 
the M{a) data, and therefore include a random noise component. The recovered spectra, 
B and S 30y g, are phase-corrected, but must be wavenumber- 2 and intensity-calibrated 2,5 (see 
also Section 5). 

From typical balloon altitudes the space spectrum <S 30 consists of a few narrow emission 
lines with no continuum emission, the limb spectra Sg contain many lines superposed on 
a substantial continuum, and the blackbody spectrum B consists of continuum alone. We 
assume that the background spectrum E(a)e l ^ a ^ is constant for the duration of an obser- 
vation sequence, and that e(cr) is constant for any given balloon flight. As will be shown in 
the next section, e does indeed seem to be constant with respect to time. 

4. Estimate of non-linear phase 

As mentioned in Section 2, each interferogram consists of a low-resolution two-sided part 
(maximum OPD range of ±L\ = ±1.2 cm) and a high-resolution one-sided part (maximum 
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OPD of L 2 = +120 cm). Both the low and high resolution interferograms are used in the 
steps which follow; we will use m to represent the transform of the low-resolution two-sided 
interferogram, and M to represent the high-resolution transform. The basic procedure is to 
estimate £30 using the high-resolution spectrum M 30 , estimate x 90 using the low-resolution 
spectrum m 9 0, and then estimate e(a) by calculating the phase of the difference m 90 — m 30 
at low resolution. In the following step-by-step procedure we have omitted some details 
(tapering functions used for low- and high-resolution transforms, zero-padding, and the 
like) for brevity and clarity; further details regarding our routine implementation of the 
fast- Fourier transform are presented in Ref. 2. 

Our method differs from one previously described 6 in that it requires two auxiliary spectra 
but works for almost totally one-sided interferograms. The radiometric calibration method 
described in Ref. 5 calibrates and phase-corrects spectra in a single step, but requires saving 
the complex transforms of hot and cold calibration spectra. Neither method 5,6 is capable 
of determining the phase factor c hu 7 $x when the two-sided transform is dominated by the 
complex background E(a). 

Step 1. To begin with, we examine the nominal (theoretical or experimental) spectra to 
be analyzed, and select a large number (40 to 60) of spectral lines which we anticipate should 
be intrinsically symmetric with respect to reflection about their central wavenumbers. For 
each line in this group, estimate a value of x 30 by symmetrizing the line in the full resolution 
30° spectrum. As can be seen by examining Eq. 1, we can determine the phase (modulo ir) 
for a spectral line at o l by finding the phase $ for which M 30 e _l/?< is symmetric near the line 
center. We calculate x 30 by fitting the function (3(a) = 27tax 30 to the derived set of (cr*, pi) 
values. The residuals from the fit are a combination of the error in estimating Pi and the 
non-linear phase term e(ai). In practice, the uncertainty in Pi is too large to get a good 
estimate of e(crj) from the residuals. 

Step 2. Calculate m' 30 , the low- resolution space spectrum with the linear phase removed: 
m 30 (a)e- i2 ™ i3 ° = [S 30 + J5e*]e* (5) 
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= m 30- (6) 

Step 3. Estimate x 90 by calculating the phase from the transform of the short two-sided 
part of the blackbody interferogram. We calculate the phase 0 t at each point in regions 
where B(a) » E(a), and again we fit a linear function to the set of (a l: Si) to derive x 90 . 
The residuals from this fit would give a good estimate of c(a), but only in regions where the 
background is small. 

Step 4 ■ Calculate m' 90 , the blackbody spectrum with the linear phase removed: 


m 90 (a)e~ i2ncri90 = [B + Ee l *}e u , 

(7) 

= ro 90 - 

(8) 

Step 5. Calculate the non-linear phase e(a): 


m' 90 ~ rn'sQ = (B + Ee^e* - (S 30 + Ee^)e { \ 

(9) 

= (B- 5 30 )e i£ ; 

(10) 

e(a) = tan' 1 “""HI 

^ 5ft ( 777-90 ~ m 3o). 

(11) 


Step 6. Repeat the process for a number of pairs of M 90 and M 30 , average the results, 
and smooth the average to produce the final estimated non-linear phase term e(a). 


In Fig. 2 we show e(cr) for two sets of measurements made several years apart using 
different beamsplitters. The estimated phase is unreliable near minima in the beamsplitter 
efficiency function because in these regions B < E and the difference m' 9Q - m 30 is dominated 
by noise. Outside these regions the two phase estimates are in excellent agreement, indicating 
that e(cr) is constant in time and independent of beamsplitter thickness. A discussion of the 
possible physical origin of the non-linear term c(cr) is found in Ref. 5. 

5. Background estimate and normalization 

Step 7. Estimate the background Ee l4> by clipping lines from the high-resolution phase- 
corrected spectrum M zo e~ li and smoothing the result, as follows. First remove the non-linear 
phase e(a) from the high- resolution space spectrum M' zo : 
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= \(S X + Ee'^e-' 1 , 
= Sjo + Ee * 


( 12 ) 

(13) 


Then effectively subtract £30 by clipping the sharp spectral lines from the result: 

- Szo = ( 14 ) 

= M" . (15) 

Clipping lines is equivalent to subtracting S 30 since the continuum emission in the space 
view is negligible. The clipping is accomplished numerically by replacing each apparent 
sharp spectral line (above a preset noise threshold) by a straight-line segment which bridges 
the adjacent spectral regions. 

Step 8. Calculate the estimated limb spectrum Sg(a) for the remaining elevation angles 
by estimating x e as for x 30 , phase correcting, and subtracting the background as follows: 

Mge~ il e~ i2waie - M" 0 = (S 9 + Ee*)^™ 19 ^ e ~* - Ee i<t> , (16) 

= Sq(g). (17) 

Step 9. Calculate the observed response B(a) to the blackbody source: 


M^e-* - M" = B(a). 


(18) 


Step 10. Calculate the intensity-calibrated limb spectra Sg al (a) from the derived Sg(a) 
and B(a): 


_ 5g(cr)P(T 66 ,cr) 

_ £(cr)P(277 K,a)’ 


(19) 


where P(T , a) is the Planck function, T bb is the measured temperature of the reference black- 
body source, and 277 K is a standard reference temperature, chosen for convention in our 
data analysis procedure. The resulting spectra S^(a) are dimensionless and independent of 
the instrument spectral response, which is convenient for plotting and visualization. Spectra 
in physical units (erg cm -2 s -1 st -1 ) can easily be recovered by multiplying by P( 277 K, a). 
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6. Results 


We judge the success or failure of the phase determination by comparing measured and 
calculated emission line spectra and looking for systematic differences. As mentioned in 
Section 4, the scatter in the estimated phase $ is too large to be able to judge the effect of 
the non-linear phase term on single spectra. In coadded spectra, however, the improvement 
is considerable in regions where e(cr) is large, particularly in the region 560-610 cm -1 where 
the non-linear term is as large as 15°. This region is used for fitting CO 2 and N 2 0 in FIRS- 
2 spectra, and correcting for the non-linear phase term eliminates a systematic difference 
between observed and calculated line shapes in this region. Spectra which have been phase 
corrected using the method outlined in this paper have emission line profiles which are visibly 
more symmetric than results from assuming a linear phase function, as shown in Fig. 3, and 
are a much better match to calculated spectra. Failing to correct for the non-linear phase 
term not only results in a shift in the line position, but also produces highly asymmetric 
line wings. 

7. Conclusion 

We describe a method for phase correcting and calibrating emission spectra from mostly 
one-sided interferograms. The method, developed for the FIRS-2, is effective in the presence 
of an out-of-phase background. The algorithm is in routine use to reduce FIRS-2 spectra, 
and has been used for several years now with good success. 
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FIGURES 



Fig. 1. FIRS-2 spectra recorded at a balloon altitude of 36 km. Vertical scale (arbitrary units) is 
the same for each graph. The elevation angles are indicated; at 90° the telescope views a reference 
blackbody. 
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Fig. 2. Non-linear phase term l estimated for two different Mylar beamsplitters. The data 
from 75-220 are from the far-infrared detector, and data from 330-700 are from the mid-infrared 
detector. 



Fig. 3. A portion of a flight spectrum which shows the effect of ignoring the non-linear phase 
term e. The upper curve was produced using the phase-recovery method outlined in this paper; 
the heavy curve indicates a calculated spectrum for the same region. The lower curve shows the 
result after assuming that e = 0. The baselines for the upper curve and calculated spectrum have 
been shifted by 0.12 for clarity. 
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APPENDIX B 


Ozone Production and Loss Rate Measurements in the Middle Stratosphere 



